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Electronic states of single-component molecular metals M(tmdt)2 (M=Ni, Au) are stud- 
ied theoretically. We construct an effective three-band Hubbard model for each material by 
numerical fitting to first-principles band calculations, while referring to molecular orbital cal- 
culations for the isolated molecules. The model consists of two kinds of base orbital for each 
molecule with hybridization between them, i.e., a 7r-character orbital for each of the two tmdt 
ligands, and, a pc/7r-orbital for M=Ni or a pdcr- orbital for M=Au centered on the metal site; 
this indicates that these materials can be considered as novel multiband it-d systems. We find 
that both orbitals contribute to realize the metallic character in Ni(tmdt)2. The origin of the 
antiferromagnetic transition observed in Au(tmdt)2 is also discussed based on this model. 
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The success in synthesizing single-component molecu- 
lar metals (SCMM) by the Kobayashis and co-workers 1 3 
opened a new route in realizing molecule-based elec- 
tronic conductors. These compounds intrinsically pos- 
sess conducting nature by self-assembly of same kind of 
neutral molecules, in sharp contrast with conventional 
molecular conductors composed of two or more chemical 
species with charge transfer (CT) between them, which 
was the unique way in realizing charge carriers in molec- 
ular crystals for a long time. The first SCMM synthe- 
sized was Ni(tmdt)2, 1 whose crystal structure is shown in 
Figs. 1(a) and (b); its electrical resistivity p(T) decreases 
as decreasing temperature down to 0.6 K and the exis- 
tence of a Fermi surface (FS) was confirmed by de Haas- 
van Alphen oscillation. 4 Among several SCMM reported 
since then, 2,3 Au(tmdt)2, 5 which is isostructural to its Ni 
analog and likewise metallic down to low temperatures, 6 
has attracted interest because of its high antiferromag- 
netic (AF) transition temperature among molecular con- 
ductors at Taf=H0 K. 5 ' 7 ' 8 In this study we focus on the 
electronic properties of these two compounds, M (tmdt) 2 
with M=Ni and Au, to shed light on the nature of the 
electronic states in this new class of materials. 

All SCMM known to date are composed of transition 
metal complex molecules with extended-TTF dithiolene 
ligands of the form M{V)2 (M: metal, L: ligand), re- 
flecting the strategy of "frontier molecular orbital (MO) 
designing" . 9 They are chosen to satisfy the desired con- 
ditions for realizing SCMM: (i) small energy separations 
AE between frontier MO in the isolated molecules, and 
(ii) existence of large inter-molecular transfer integrals 
when crystals are formed. As discussed in ref. 9 in de- 
tail, M(V)2 molecules are suitable for condition (i), since 
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Fig. 1. Crystal structure of M (tmdt) 2 (M=Ni, Au) [(a) and (b)], 
and a schematic view of the model discussed in this paper where 
the lattice sites are represented by circles [(c) and (d)]. In (b) 
and (d), some of the intermolecular and intersite bonds with large 
transfer integrals (see text) are shown. In (c), two-dimensional 
plane of ligands projected onto the 6c-plane is represented in 
light blue, whose network within the a6-plane is shown in (d). 

their frontier MO are usually characterized by bonding 
and anti-bonding combinations of "virtual" frontier MO 
of the two ligands. Then AE is roughly twice the effective 
transfer integral between them, which can be tuned small 
by using appropriate L and M. As for condition (ii), 
TTF-type skeletons produce large inter-molecular over- 
laps efficiently as heavily used in CT salts. When these 
conditions are fulfilled for molecules with even number 
of electrons, widths of the bands near the Fermi energy, 
ep, originated from frontier MO can exceed AE, leading 
to band overlaps so that FS would exist, i.e., resulting in 
semi- metallic states. As for molecules with odd number 
of electrons (which are rather rare) the latter condition 
is necessary to avoid localization due to electron correla- 
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Fig. 2. Band structures of (a) Ni(tmdt)2 and (b) Au(tmdt)2 
near the Fermi energies, between r(0,0,0)-point, and X(7r,0,0), 
Y(0,tt,0), Z(0,0,tt), S(7r,7r,0), U(7r,0,7r), T(0,7r,7r), and R(7r,7r,7r)- 
points (black lines). Fitted results and the Fermi levels are shown 
by red solid and blue dotted lines. 

tion, typically Mott insulating states. 

The band structures of Ni(tmdt)2 and Au(tmdt)2 have 
been calculated by the extended Hiickel tight-binding 
scheme 1 as well as by ab-initio first-principles calcula- 
tions. 10 12 The results show that several bands originated 
from different MO are situated near ep with apprecia- 
ble mixings between them, consistent with the above- 
mentioned strategy. Note that this is very different from 
most of the CT-type molecular conductors where the 
bands near ep are formed by one 7r-character frontier 
MO per molecule. There, effective Hubbard-type models 
have been successfully applied to study their low-energy 
physical properties theoretically. 13 The purpose of this 
work is to investigate whether we can make use of such 
effective models for SCMM as well. In the following, by 
comparing the first-principles band structures and MO of 
isolated molecules, we construct a three-band Hubbard 
model based on virtual MO of the molecules, elaborating 
the above-mentioned two ligands MO picture. 9 Further- 
more, we investigate this model to pursue the origin of 
the AF transition observed in Au(tmdt)2, which has the- 
oretically been proposed by Ishibashi et al. 11 to be due 
to FS nesting based on first-principles calculations. 

These compounds are nominally written as 
Ni 2+ (tmdt-) 2 and Au 3+ (tmdt 15 -) 2 , with metal 
ions having (3d) 8 and (5d) 8 configurations. 2 ' 3 Since 
neutral tmdt has even number of electrons, the Ni (Au) 
salt has even (odd) number of electrons per molecule. 
In Fig. 2, their band structures near ep calculated by 
first-principles method are shown. 11 ' 12 Due to the re- 
semblances in their electronic configuration and crystal 
structure, overall features of the four bands shown in 
the figure are similar between the two compounds. A 
difference is seen in the top band which is isolated in 
Ni(tmdt)2, whereas it is close to the next band near the 
X and S points in Au(tmdt)2. In contrast, in the case of 
Ni(tmdt) 2 the bottom band almost touches to the third 
band from the top in several regions of the Brillouin 
zone, while in Au(tmdt) 2 it is separated from the others. 
Therefore three bands each are relevant near ep, i.e., the 
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Fig. 3. Molecular orbitals from which the bands near the Fermi 
energy are formed. 12 The rough spatial extensions of the L and 
M orbitals explained in the text are indicated. 

lower (upper) three bands for the Ni (Au) salt, which we 
will consider in the following. The even/odd electrons 
per molecule in the two compounds result in electron 
and hole pockets with equal volume in Ni(tmdt)2 and a 
half- filled band in Au(tmdt)2. 

These bands near ep originate from three MO for each 
molecule, which are shown in Fig. 3. 12 A close resem- 
blance is seen between HOMO of Ni(tmdt) 2 and SOMO- 
1 of Au(tmdt) 2 which, as discussed in ref. 9, have the 
bonding nature between the two ligands MO. On the 
other hand, in the corresponding anti-bonding orbitals, 
i.e., LUMO of Ni(tmdt) 2 and SOMO of Au(tmdt) 2 , the 
contributions from the d-orbital at the metal site and 
the p-orbital at the four surrounding S atoms are more 
prominent in the former than in the latter. 12 Consider- 
ing these characteristics, we attempt to reconstruct MO 
in Fig. 3 by three parts of the molecules, as follows. In 
Ni(tmdt) 2 , we can postulate that LUMO, HOMO, and 
HOMO-1 are approximately described by anti-bonding, 
non-bonding, and bonding combinations of three vir- 
tual MO within the molecule, 0li — </>L2 + Cab0M(pd7r)> 
</>li +</>L2, and L i -0L2 -Cb^M(pdTr), respectively. Here, 
0li and (/>L2 are the two "ligand (L) orbitals" close to 
LUMO of the embedded TTF molecule, and </>M(pd7r) is 
the "metal (M) orbital" centered at the Ni site, and 
c a b and Cb are some coefficients. Note that we consider 
this M orbital not as a bare 3d orbital of Ni as consid- 
ered in ref. 9, but rather close to the anti-bonding pdir 
orbital between the 3^^-orbital 14 of the Ni atom and 
the four 2p-orbitals of the surrounding S atoms. As for 
Au(tmdt) 2 , LUMO is another "M orbital" roughly be- 
ing the anti-bonding pda orbital, 0M(pd<r)> between the 
Au 5^^-orbital 14 and the surrounding S 2p- orbitals. 12 
On the other hand, by neglecting the contribution from 
0M( P dn) smaller than in Ni(tmdt) 2 , SOMO and SOMO-l 
of Au(tmdt)2 are described as 0li — </>L2 and 0li + </>L2- 
We stress here that the occupation probability of the 
d-orbitals itself is small in both M orbitals, but their 
characters show up through hybridization with the S 2p- 
orbitals. 

Then, assuming the above two L and one M orbitals, 
namely, Ll, L2, and M^pdiv) for Ni(tmdt) 2 or M(pda) for 
Au(tmdt) 2 , as base orbitals, we can construct an effective 
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Table I. Fitted values of tij and A in eq. (1) with Uj j =Um = 
U = 0.3 eV. For example, "A[100] Ll-M" denotes Uj between 
LI orbital belonging to M(tmdt)2 at the origin and M orbital of 
M(tmdt) 2 at (a, 0, 0). 





site pair 


Ni(tmdt)2 


Au(tmdt)2 


intra-mol. 


L1-L2 


11.3 meV 


-54.3 meV 


A[100] 


Ll-Ll, L2-L2 


-80.8 


-95.5 


A[100] 


L1-L2 


-7.8 


-5.4 


B[lll] 


L1-L2 


216.0 


208.0 


0[1U1J 


L1-L2 


9.7 


13.0 


Q[001] 


T1T1 TOTO 

L/l-L/1, LiZ-LiZ 


9.0 


3.6 


Q[001] 


T 1 TO 


123.0 


108.7 


P[211] 


T 1 TO 

Ll-Lz 


39.8 


37.4 


K|U11J 


T 1 TO 


9.7 


A O 

4.Z 


intra-mol. 


T 1 A /T l\TTO 

L/l-M, NL-LiZ 


-218.0 




A[100] 


Ll-M, M-L2 


14.9 


38.3 


A [100] 


M-Ll, L2-M 


14.3 


1.0 


B[lll] 


Ll-M, M-L2 


41.9 


28.0 


C[101] 


Ll-M, M-L2 


-21.1 


-17.9 


Q[001] 


Ll-M, M-L2 


26.1 


11.6 


R[011] 


Ll-M, M-L2 


1.8 


17.2 


A[100] 


M-M 


-21.5 


94.3 


B[lll] 


M-M 


22.3 


-1.0 


C[101] 


M-M 


-0.6 


1.8 


Q[001] 


M-M 


-3.3 


1.5 


A 




11.5 


717.7 



three-band Hubbard model for each material written as 
H = Uj (4s c js + h-c.) + A Y, («it + nn) 

ijs i(M) 
i(L) i(M) 

where c\ s is the creation operator of an electron with spin 
s(=T/l) at the ^h site, and n^ s = cj s c is . Uj are trans- 
fer integrals between site pairs i and j, and [/"l and /7m 
are on-site Coulomb energies for the L and M orbit als, 
respectively. The one-particle energy difference between 
the two kinds of orbital is denoted as A, which is mea- 
sured from the L orbital. The summation Y]-- e is taken 
for all site pairs considered (L-L, L-M, and M-M pairs), 
while J2i(L) an d J2i(M) are f° r s ^ es belonging to the L 
and M orbitals, respectively. 

To find out the appropriate parameters in eq. (1), 
we perform numerical fitting of the energy dispersion 
for its paramagnetic Hartree-Fock solution, in order to 
maintain consistency with the spin-dependent calcula- 
tions below, to the first-principle band structures. We 
set Ul=Um = U for simplicity, which are considered to 
be of the same order deduced from similar spatial ex- 
tensions of the two orbitals seen in Fig. 1. The fitted 
results for U = 0.3 eV are shown in Fig. 2, which show 
good agreement with the first-principles band structure, 
particularly near ep. 

The fitted parameters are listed in Table I. Uj for the 
L-L pairs are all similar between the two compounds, 
naturally expected since they are isostructural. There 
are several inter-molecular Uj with absolute values much 
larger than that for the intra- molecular L1-L2 pair, which 
form a two-dimensional network in the (001) plane as 
shown in Figs. 1(c) and (d). There is a large Uj for the 
bond along B[lll] between the different ligands (L1-L2), 
of about 0.2 eV which is twice larger than the others; 
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they form dimers. 

The transfer integrals for L-M and M-M pairs show dif- 
ferent aspects between the two compounds. In Ni(tmdt)2, 
the intra-molecular L-M pair with \Uj \ ~ 0.22 eV, which 
is comparable to the above B[lll] L1-L2 bond, brigdes 
the two-dimensional L networks, and make the elec- 
tronic structure three-dimensional. On the other hand, 
for Au(tmdt)2, the intra-molecular L-M bond is not in- 
cluded in the fitting due to the tt and a symmetry of 
the two orbitals; Uj for inter- molecular L-M pairs are 
rather small but responsible for the three-dimensional 
ir-pda hybridization. 

The L-M energy difference A in Ni(tmdt)2 is surpris- 
ingly small, i.e., the two kinds of orbitals are almost de- 
generate in energy, indicating that the three bands are 
formed by completely mixed L and M orbitals. This sug- 
gests that both orbitals play crucial roles in forming the 
bands crossing ep, and therefore, in realizing the metallic 
state in this compound, in contrast with the two ligand 
MO picture as has been discussed in the litteratures. 9 ' 10 
As for Au(tmdt) 2 , the fitted value for A = 0.72 eV sug- 
gests that the top band is formed mainly by the M(pda) 
orbital. 11 We note that the difference in the value of 
U used in the fitting barely affects the parameters for 
Ni(tmdt)2, but for Au(tmdt) 2 A is affected while Uj are 
not; A should be considered as an effective L-M energy 
difference. 

Next we will proceed to the analysis of the model above 
for Au(tmdt)2 and discuss its magnetic instability. In 
the following we investigate the electronic properties of 
this model by varying A from its fitted value (= Ao) 
while fixing other parameters, i.e., Uj as in Table I and 
[7=0.3 eV. Within spin-dependent Hartree-Fock approx- 
imation, we have sought for magnetically ordered solu- 
tions with commensurate AF spin patterns up to 2x2x2 
of the original unit cell, and a stable solution consistent 
with that obtained in ref. 11 is found. This state has 
anti-parallel magnetic moments on the two L orbitals 
within each molecule, and characterized by the wave vec- 
tor q = (a*/2,0,0), whereas no moment appears on the 
M sites (see Fig. 8 in ref. 11 and the inset of Fig. 4 (a)) In 
Fig. 4(a), the calculated transition temperature Taf for 
this solution as well as its magnetic moment at T = 0, 
mo, are shown as a function of A. The magnetically or- 
dered state has small but finite density of state at ep, 
D(ep), as seen in the inset of Fig. 4 (c), i.e., the sys- 
tem is still metallic. 11 Taf and mo show a maximum at 
around A = Ao, whose origin is understood as follows. 

In Fig. 4 (b) and (c), we show the FS and the density 
of states D(e) of the paramagnetic solution, for A = Ao 
and Ao ± 0.15 eV. One can see that FS changes little 
when A is increased from Ao, and its nesting property 
found in ref. 11 is maintained; therefore the nesting is due 
to the 7r-network of L orbital. When decreasing A from 
large values, the slight increase of Taf and mo at A ^ Ao 
is due to the enhancement of D(ep) by the increase of L- 
M hybridization as seen in Fig. 4 (c). This gives rise to 
an increase of the whole profile of the bare susceptibility 
Xo(q) (not shown; see Fig. 6 in ref. 11), directly related 
to the magnetic instability. Further decrease of A actu- 
ally raises D(ep) more rapidly, but leads to modification 
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Fig. 4. (a) Transition temperature Taf and T = magnetic mo- 
ment mo, for the antiferromagnetic (AF) state, (b) Fermi surface 
(the Brillouin zone is deformed to a rectangular solid), and (c) 
density of states -D(e), for the paramagnetic metallic state in the 
Hartree-Fock calculation for Au(tmdt)2 by varying A (the inset 
shows change in D(e) near ep between the paramagnetic and AF 
states for A = Ao). A schematic view of the AF ordered spin 
pattern of the L orbital sites in the a6-plane is shown in the inset 
of (a), whereas M sites have no moment. 

of FS and easily destroys the nesting, which results in 
the sudden decrease in Taf and mo- We note that we 
find a shift of peak position in Xo(q) to an incommen- 
surate wave vector q = (a*/2, 5, 0) in a narrow region of 
0.59 eV < A < 0.63 eV, then an incommensurate AF 
could be stabilized there which is not considered in our 
calculations. Still we expect that the rapid decrease of 
Taf holds even when such possibility is included. 

Let us compare our results with the experiments on 
Au(tmdt)2. If the transition at 110 K seen in the experi- 
ments is due to the spin-density- wave (SDW) formation 
by FS nesting, a signature in p(T) should appear due 
to the change in the band structure, as a hump in p{T) 
below the transition temperature as found in metallic 
Q r 16 This was no t seen in previous measurements on 
compressed pellet samples 5 nor on microcrystals, 6 while 
a small kink has recently been detected by a prelimi- 
nary measurement on a single crystal 15 ; further mea- 
surements on single crystal samples are awaited. A pos- 
sible origin of the high Taf among molecular conductors 
is the three-dimensionality of the system as has been dis- 
cussed in ref. 8. This is not realized in conventional CT 
salts showing SDW, e.g., in TMTSF 2 X, due to the exis- 
tence of the anions or cations providing carriers into the 
7r network resulting in low-dimensional electronic struc- 
tures, which lower transition temperatures in general due 
to fluctuation effects. Our spin-dependent Hartree-Fock 
results, which can be interpreted as such SDW state, ac- 
tually show the same order of Taf to the experiments. 

However, as for the magnetic moments in the ordered 
state, analysis of a recent ^-NMR experiment suggested 



that they are rather large, more than 0.6 pb per tmdt, in- 
consistent with our results showing saturated moments 
at T = less than 0.2 p^ per tmdt. Such large mag- 
netic moments in molecular conductors are reminiscent 
of strongly correlated systems such as Mott insulating 
states in the K-ET2X family 18 or charge ordered states 
in the TMTTF 2 X and <9-ET 2 X families. 19 In fact, the 
AF ordered spin pattern we disscussed in this paper as 
shown in the inset of Fig. 4(a) can also be considered as 
"AF between dimers" , 20 i.e., ferromagnetic within dimers 
formed by the B[lll] L1-L2 pairs and AF between them, 
and an insulating gap opens when larger U ( ^ 0.4 eV) is 
used in the calculation. This state can be viewed as the 
dimer-Mott insulator, 13,20 although it cannot be identi- 
fied within Hartree-Fock approximation. In the above- 
mentioned strongly correlated states, experimentally the 
systems show insulating behavior at low temperatures 
and p(T) shows no anomaly at the magnetic transition 
temperatures; these are both distinct from the experi- 
mental data for Au(tmdt) 2 with metallic behavior and 
the kink at Taf in p(T). Nevertheless, the high Taf and 
large magnetic moments seen in the experiments can be 
due to the closeness to such dimer-Mott insulating states, 
whose investigations are left for future studies. Finally, 
we note that a low transition temperature below 10 K 
was recently found in Au(tmstfdt) 2 , 17 a selenide analog 
of Au(tmdt) 2 , despite they are isostructural and have 
similar electronic structures. This might be a result of 
the fragileness of the FS nesting by L-M hybridization, 12 
seen in our calculation. 

In summary, we have constructed an effective three- 
band Hubbard model for single-component molecular 
metals M (tmdt) 2 (M=Ni, Au) which indicates that this 
system is a multiband ir-pdn / pda system. We have found 
that the close degeneracy of 7r and pdn "virtual" MO 
states in Ni(tmdt) 2 is a key to realizing its metallic state. 
The origin of the antiferromagnetic transition observed in 
Au(tmdt) 2 is discussed based on spin-dependent Hartree- 
Fock calculations. We believe that this new category of 
materials will be promising for finding richer phenomena 
in molecular crystals. 
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